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We study a chiral relativistic theory of nuclear matter aimed at the desciption of both the binding 
and saturation properties and the QCD properties, quark condensate and QCD susceptibilities. For 
this purpose the nucleon scalar response of the quark-meson coupling model is introduced in the 
linear sigma model. The consequences for the nuclear and the QCD scalar susceptibilities are 
discussed. 
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I. INTRODUCTION 

The order parameter associated with spontaneous breaking of chiral symmetry, namely the chiral quark condensate, 
is influenced in the nuclear medium by the mean scalar field which possesses the same quantum numbers 1]. Its 
fluctuations in the medium as well are intimately related to the in-medium propagation of the scalar field On 
the other hand the properties of the scalar field influence in a crucial way the question of the nuclear binding. It 
is therefore important for both problems to reach a description of the nuclear dynamics which satisfies the chiral 
constraints and is able to correctly reproduce the binding and saturation properties. The key ingredients can be 
found in the theory of quantum hadrodynamics (QHD) which, however, in its original form lacks the chiral 

aspects. We have recently shown that, starting from a chiral effective model, it is possible to give a chiral status 
to QHD by imposing a chiral invariant character to the scalar field of QHD 1|. However there is a well identified 
problem concerning the nuclear saturation with usual chiral effective theories 5]. Independently of the particular 
chiral model, in the nuclear medium one moves away from the minimum of the vacuum effective potential (mexican 
hat potential), i.e., into a region of smaller curvature. This single effect, equivalent to the lowering of the sigma mass, 
destroys the stability, creating problems for the applicability of such effective theories in the nuclear context. The 
problem can be cured with the introduction of the nucleonic response to the scalar field which is the central ingredient 
of the quark-meson coupling model, QMC 0,01- Indeed, the nucleon polarization under the infiuence of a scalar field 
implies a density dependence of the scalar nucleon coupling constant, which actually decreases with increasing density. 
This effect can counterbalance the decrease of the curvature of the mexican hat potential and restore saturation. The 
QMC model is quite successful in the phenomenology of nuclear matter or finite nuclei. It is the aim of this work 
to incorporate its main concept, i.e., that of a nucleonic response to a scalar field, in an effective chiral theory, the 
linear sigma model. This concept is the only feature that we retain of QMC, which goes beyond since it evaluates the 
scalar polarizability kjvs in a bag model. For us kjvs is a parameter that we can adjust. We are however guided by 
the values obtained in the original model. In particular the sign found in QMC is crucial. The fact that the nucleonic 
scalar response manifests itself at the nuclear level is plausible. This is familiar in other situations. For instance under 
the influence of an isovector axial field, the nucleon converts into a Delta. This axial polarizability of the nucleon 
has several physical implications (such as the in-medium renormalization of the axial coupling constant). What is 
unusual in the scalar polarizability is its sign. In QMC, the scalar polarizability originates from relativistic effects of 
quarks confined in a bag (Z graphs) and it is positive. It screens the scalar field (diamagnetic effect). 

Inside of our framework, i.e., the linear sigma model with a nucleonic scalar response, we will choose the parameters 
so as to be compatible with the known saturation properties of nuclear matter: saturation density, binding energy 
and compression modulus. The resulting density dependence of the sigma mass becomes quite mild while, in the 
absence of the nucleonic scalar response, this mass drops rapidly with increasing density, preventing saturation to 
occur. We also derive the QCD scalar susceptibility of the nuclear medium. An interesting by-product of our work 
is the obtention of the residual particle-hole interaction in the isoscalar channel. It arises from in-medium sigma 
and omega exchanges. The introduction of the nucleonic scalar response leads to a rapid density dependence of this 
interaction. It turns from attraction to repulsion near saturation density. Accordingly the collective nature of the 
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response undergoes a qualitative change, from a softening effect at low density to a hardening one above saturation 
density. 

Section 2 is devoted to the formal derivation from the equation of state of various quantities: sigma mass, QCD 
susceptibility, nuclear response. In section 3 we perform the numerical evaluations of these quanties, ensuring the 
compatibility with the nuclear phenomenology and we discuss the consequences. 



II. EQUATION OF STATE AND SCALAR FLUCTUATIONS 

We start with the Lagrangian introduced in our previous work I'l , dropping its pionic part which does not contribute 
at the mean field level: 

L = ia^sa^s - V{s) + iN-f^'^^N - Mn (1 + ^ ) NN (1) 



where s is the chiral invariant field associated with the radius 5 = s + of the chiral circle. V{s) = Vo(s) — cS is the 
vacuum potential which can be split into Vq{s) = (A/4)((/7r + s)^ — v^^ responsible for spontaneous chiral symmetry 
breaking and the explicit symmetry breaking piec, —cS, where c = /ttW^- As usual in QHD we add a coupling to an 
omega field : 

Cv = Idf^ajd^LO + ]^mlu^ ~ g^Lo N^N. (2) 
As for the nucleonic response, kats, to the scalar field s, we incorporate it in the following extra term in the lagrangian : 

C^^-^nNss" NN. (3) 

The susceptibility kjsis embeds the influence of the internal nucleon structure. We will discuss later its possible s 
dependence. At the mean field level the energy density is given by : 

^ Q{pp - p) E;{-s) + V{s) + ^p^ (4) 

where the baryonic density is related to the Fermi momentum through : 



and E*{s) — y^p-^ + M'^{s) is the energy of an effective nucleon with the effective mass : 

= Mn (l + ^) + \ '^NS s'. (6) 

The expectation value, 5 = /jr + s, of the S field, plays the role of the chiral order parameter. It is obtained by 
minimizing the energy density : 

^_=g*ps + V'is)=0 (7) 
with the following expressions for the scalar density, ps, and the scalar coupling constant g^ : 

Notice that the density dependence of g^ entirely arises from the susceptibility term. Since the mean scalar field is 
negative and the sign of kns positive, gg is a decreasing function of the density. In the vacuum the scalar coupling 
constant of the model is gs = M^/ /tt- The in- medium sigma mass is obtained as the second derivative of the energy 
density with respect to the order parameter : 



2N 

U ' 2 V./J / ' '^"""^ ' ds 



"4 l + ^ + ^f^) ]+-NsPs+9*s^ (8) 



where in the second Hne we have taken V" (s) in the chiral Umit. The mean scalar field s being negative, the term linear 
in s (which appears from the curvature of the effective potential) in itself lowers the sigma mass by an appreciable 
amount (~ 30 % at po). This is the chiral dropping associated with chiral restoration emphasized by Hatsuda et al 
01 , which they suggested to be the origin of the strong medium effects found in 2-k production experiments "To', Till . 
However the scalar susceptibility term (second term in kjvs on the r.h.s of eq. ISJ which is numerically important, 
counterbalances this in- medium mass dropping. Indeed in QMC, the sign of kns is positive, an essential element of 
the model for reaching saturation. The influence of the nucleonic scalar susceptibility thus increases the sigma mass, 
counterbalancing the chiral restoration effect. In pure quantum-hadrodynamics instead the chiral softening of the 
sigma mass is ignored and saturation is obtained through a large polarisation of the Dirac sea of nucleons. Hence the 
sigma mass evolution has to be discussed in connection with saturation properties. As for the last term {g*g dps / ds) , 
it also writes : 

* dp*s *2 [ 4rf^p „ , . 
^^^-3^ J (2^®(^^-P)£f 

We have shown in a previous work 12] that it actually corresponds to the nuclear response associated with NN 
excitation. In practice it is small and it can be omitted. 

We will now derive the in-medium chiral condensate and the QCD scalar susceptibility. They are related to the 
first and second derivatives of the grand potential with respect to the quark mass m at constant chemical potential 
p. The baryonic chemical potential is obtained as : 

de 0^ I 

p= — ^E*p^^p with E*p = Jpl + M*^(-s) (10) 
op V 

from which one deduces that the baryonic density is controled by the chemical potential according to : 

/. = /^e(M-i.;-|-.) (11) 

while the scalar density writes : 
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The grand potential, which is obtained through a Legendre transform, can be written in the following form: 



uj{p) = e - pp 
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('^; + f - ") - i=; - 1;.) + ViA - (13) 

Notice that the minimization (equation [TJ can be equivalently obtained from the condition [duj/ds)^ = 0. In order 
to derive the condensate and susceptibility we point out that, in the context of this model, what plays the role of 
the chiral symmetry breaking parameter is the quantity c = /^mj which enters the symmetry breaking piece of the 
potential. Hence : 

where we hase used the Feynman-Hellman theorem and the explicit expression of dc/dm given by the model to leading 
order in the quark mass to. Accordingly the in-medium scalar susceptibility is given by : 



The derivative (dS/dc)^ is obtained by taking the derivative of the minimization equation Q with respect to the 
parameter c. This gives : 



with no(0) = —2M'^pp/'k^. Notice that no(0) is nothing but the non-relativistic free Fermi gas particle- hole polar- 
ization propagator in the Hartree scheme, at zero energy in the limit of vanishing momentum. The derivative of the 
baryonic density is obtained by taking the derivative with respect to c of eq. 1111 with the result : 
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It follows that : 



ds\ Lnss(o) — 

dc J ^ m*^ m*^ m*^ 



(18) 



where 1155(0) is the full scalar polarization propagator (in which we include the coupling constant) 
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We now comment these results. Firstly notice that the quantity {dS/dc)f^, as written in eq. ^1 is the in-medium 
scalar meson propagator (up to a minus sign) dressed by NN~-^ excitations. Moreover it is satisfactory to realize that 
the expression of 1155 (eq. I19|l coincides with the one derived from the RPA equations in the ring approximation. In 
RPA the residual interaction can be due either to the scalar meson exchange or to the vector one. In the latter case 
the mixed propagator Hsv enters. The RPA equations read : 
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One recovers our result of eq. 1191 since, as shown in the appendix, it is possible to check that: 



(21) 
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(22) 



In our approach the Landau-Migdal parameter Fq enters the RPA denominator of eq. I19l which. at poj writes 1 -\- Fq. 
This equation also shows that our residual particle- hole interaction is density dependent, in particular through the 
density dependence of g*g and m* . Our approach provides a consistent relativistic frame in a chiral theory to derive 
this density dependence. A similar expression has been given in ref. 

As for the full vector polarization propagator Ilyy, solution of the RPA equations, it is given by : 
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Notice that the same RPA denominator enters both expressions of 1155 and Ilyy , with omega and sigma exchange on 
the same footing. The quantity rtyy, i.e., the response to a probe which couples to the nucleon density fluctuations, 
is related to the nuclear compressibility. We introduce the incompressibility factor K of the nuclear medium, defined 
as : 
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The minimization (equation0) establishes the dependence of s with respect to p. By taking its derivative, one gets 
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It follows that : 

UvvjO) _ _9p 



(26) 



This relation is well known in the non relativistic situation. We have shown that it also applies in the relativistic case 
with Hvv given above (eq. I23|l . At low densities where relativistic effects are small, there is no distinction between 
the scalar and vector propagators and we have Ilss{0)/9s^ — ^vv{0)/ 9^- More generally the relation is : 



9*3' \E*J 9l ■ 



(27) 



With our values of parameters the parenthesis represents a reduction of about 10% at po- The nuclear physics 
information on the nuclear matter compressibility leads to a value of K in the range 200 — 300 MeV. This is close to 
the free Fermi gas value, which is 230 MeV for = Mm- Thus the nuclear phenomenology which constraints the 
residual interaction at saturation density, also constraints scalar quantities, such as the scalar nuclear response and 
the QCD scalar susceptibility. 

III. NUMERICAL RESULTS AND DISCUSSION 
A. Sigma mass 

The first consequence we wish to discuss is the much debated problem of the density dependence of the scalar 
meson mass. Its explicit form (in the chiral limit) is given by eq. |S1 and we now come to the quantitative estimate. 
It is possible to get a first evaluation of the nucleon response, kjsis, from the parameters of the QMC model, as 
given in (21 . It gives for the dimensionless parameter C — {f^/2MN)KNS the value C — 0.45. Two independant 
parameters remain to be fixed : 9slTna and 9u]/mi^. Since gs = Mtv/Ztt is given by our model and since we take for 
the omega mass the vacuum value = 783 MeV, the parameters to be fixed are in fact mg- and gui- In the first 
step of the discussion we keep for the scalar potential, Vb(s), only the quadratic part, namely, Vo(s) = m^s^/2, as 
in the original Walecka or QMC model which both ignore the chiral softening of the sigma mass. In such a case the 
saturation properties (po = 0.17/to~^, E/A = —15.3 MeV) are obtained by taking = 715 MeV and g^j = 7.47. 
The corresponding incompressibility is K = 260 MeV and the effective nucleon mass at po is = 760 MeV. In 
general with the introduction oi n^s the effective nucleon mass is closer to the free value than in pure Quantum 
Hadrodynamics where ~ 600 MeV. 

If instead we introduce the full chiral Vb(s), the effect of chiral restoration associated with the dropping of the 
sigma mass (see eq. |H1 and the accompanying discussion) tends to destroy the saturation mechanism. One possible 
way to recover saturation is to increase the nucleonic susceptibility kns (or equivalently the dimensionless parameter 
C). Indeed with a larger value, C — 0.8, and keeping the other parameters to the same values we find saturation at 
po- However the binding energy E/A — —17.8 MeV and the incompressibility K — 360 MeV are too large. A slight 
readjustment of the parameters C = 0.85, nia — 750 MeV and guj = 7 can solve the problem for the binding energy 
(see dot-dashed curve of figure 1) but the incompressibility remains too large (above 300 MeV). The fact that we 
do not fully reproduce with this simple model the saturation properties is not surprising. We have assumed a linear 
response of the nucleon {kns is independent of s), which is not founded. It is easy to introduce a field dependence. 
We take the simplest linear one, imposing the vanishing of the susceptibility at chiral restoration (s = ~/ir) : 

KNS KNS{S) = ^^r, = KNS 1 1 + — ) ■ (28) 

Accordingly the effective nucleon mass becomes : 



M*^ = Mn + ^ns (^1 + ^ j . (29) 

In this case the set of parameters : g^^ ~ 6.8, = 750 MeV and C = 1, leads to correct saturation properties, 
with an incompressibility value K — 270 MeV (full curve of figure 1) and M^{po) = 760 MeV. Notice that in this 
case the corresponding value of the residual interaction Fq at po is practically zero {Fq ~ 0.03). This is due to a 
delicate cancellation between two large terms : the omega exchange and the in-medium modified scalar exchange, the 
magnitude of each term being of the order of 3. 
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P/Po 

FIG. 1: Binding energy of nuclear matter. The dashed line corresponds to the set of parameters g^) — 7, ma = 750 MeV and 
C — 0.85 in the absence of the density dependence of the nucleon susceptibility. The full line corresponds to a set of parameters 
Qui = 6.8, rricr ~ 750 MeV and C = 1 with an explicit field (density) dependence of the nucleon susceptibility as explained in 
the text. 



Coming back to the sigma mass evolution, the general behaviour is to a large extent independent of the exact field 
dependence of the susceptibility. Fig. 2 represents this evolution without (dot-dashed curve) and with the effect of the 
nucleonic scalar response, introducing or not its field dependence. The two curves with this nucleonic scalar response 
are fiat. The nucleon reaction largely suppresses the strong softening due to chiral restoration which, if taken alone, 
would not be compatible with saturation properties. A similar conclusion was reached in ref. j5(|. 



B. Nuclear responses 



Our second point concerns the nuclear response to a probe which couples to the nucleon density, scalar or vector, for 
which we have given the relativistic expressions (eq. I19I23|) . At po we have seen that the value of the incompressibility 
K of nuclear matter lies in the range 200-300 MeV. This is close to the free Fermi gas value, which is K — 230 MeV 
for = Mm- In QMC where the effective mass does not differ so much from the free one, the experimental value 
of K implies a small residual force at po- Our choice of parameters takes this constraint into account. This smallness 
of Fq results from the accidental cancellation between omega and sigma exchanges. However this delicate balance 
is upset with a change of the density. Several factors are responsible for this phenomenon. Firstly the relativistic 
factor [M^/Ep) lowers the sigma exchange while it enhances the omega one, by a density dependent amount. In 
QHD, where the effective nucleon mass is small, this infiuence is large. In QMC instead, with larger mass values this 
is not the main effect. The main one is due to the action of the nucleonic reaction kjvs which is responsible for the 
decrease of the scalar coupling constant with increasing density. The sigma effectively decouples from the nucleon 
at large density, leaving the repulsive omega interaction to dominate. The reverse occurs at smaller density with 
the increase of gs- The sigma attraction then fully develops and dominates the repulsive omega exchange. Thus, 
with increasing density, the residual interaction turns from attraction into repulsion. Each component being large the 
evolution is fast. For instance, with our parameters the resulting attraction is strong enough to produce a singularity 
of the polarization propagator at a density p ~ 0.6po- We identify it with the spinodal instability. In the density 
region below the saturation one, the responses, scalar or vector, are collective with a resulting softening of the the 
response. It is refiected as en enhancement in magnitude of the nuclear susceptibility (proportional to 11(0)). Above 
Pq instead collectivity hardens the response and decreases the susceptibility. 

The density dependence of the residual interaction is well established |0, . It is small in the nuclear interior 
and strongly attractive at the nuclear periphery. On the other hand the collective nature of the vector response is 
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FIG. 2: Density evolution of the sigma mass. Dotted line : in the absence of the field (density) dependence of the nucleon 
susceptibility with values of the parameters g^j = 7, = 750 MeV and C — 0.85. Full line : with density dependence of the 
nucleon susceptibility with g^i = 6.8, rua = 750 MeV and C = 1. Dot-dashed line :it corresponds to the case where only the 
chiral softening is included, without the effect of the nucleon susceptibility. 
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FIG. 3: Density evolution of the QCD susceptibilities normalized to the vacuum value of the scalar one calculated with the 
field dependence of the nucleon susceptibility. Full curve : scalar susceptibility. Dashed curve : difference between scalar and 
pseudoscalar susceptibilities. 



8 



confirmed by data, from (e, e') scattering, on the longitudinal response of various nuclei at small momenta. The 
charge response is a sum of a vector-isoscalar response and a vector-isovector one. Only the first part is relevant for 
our discussion. The data at small momenta (200-300 MeV/c) display a strong softening effect with respect to the free 
response. Alberico et al have attributed this feature to the collective caracter, at the nuclear periphery, of the 
isoscalar part of the charge response, in accordance with the present views. The density dependence of the residual 
force with the non-trivial evolution from attraction into repulsion with increasing density naturally follows in our 
description with sigma and omega exchange and with the incorporation of the nucleonic scalar response, necessary to 
favor the omega contribution relative to the sigma one with increasing density. 



C. QCD susceptibilities 

It is normal that the collective character of the nuclear response shows up in the QCD susceptibility, as the quark 
density fluctuations are coupled to the nucleonic ones. In ref. |2| we have derived, in the linear sigma model, the 
following relation between the in-medium QCD scalar susceptibility and the vacuum one : 



XS,va 



For a quantitative evaluation we calculate the r.h.s. in the model which incorporates the nucleonic response with its 
field dependence. The corresponding parameters are those given previously, which reproduce the saturation properties. 
The resulting density evolution is shown in fig. 3. At po the enhancement over the vacuum value is a factor of about 5. 
At lower density it becomes even larger due to the collectivity of the polarization propagator Xlgs, with the attractive 
p-h force which enhances its magnitude. At larger densities instead the p-h force becomes repulsive which shows up 
in the gradual disappearance of the medium effect, as depicted in fig. 3. 

It is interesting to contrast the evolution of the scalar susceptibility with that of the other QCD susceptibility, the 
pseudoscalar one. This last quantity is linked to the fluctuations of the quark pseudoscalar-isovector density. The 
question is wether the increase in magnitude of the scalar susceptibility due to the coupling to the nucleon-hole states 
leads to a convergence effect between the two susceptibilities, which would be a signal for chiral symmetry restoration 
since the two susceptibilities become equal in the restored phase. In ref. we have demonstrated that the evolution 
of the pseudoscalar susceptibility follows that of the quark condensate, with : 

XPS = (31) 

which diverges in the chiral limit, as it should. The condensate on the r.h.s. is a function of the density and also of 
the quark mass. Since we are not considering here pion loops which would produce a non- analytical behavior in ruq, 
far from the restoration density we can expand the condensate, writing : 

imip)) {qqip^mq ^ o)) + xsm^ 

Xps = = ■ [oZ) 

rriq niq 

This equation provides in fact the difference between the pseudoscalar and the scalar susceptibilities. The important 
point is that this difference which is governed by the value of the condensate in the chiral limit is a smooth fuction of 
the density, which does not feel the large fluctuations near the spinodal instability. This means that the large increase 
of the scalar susceptibility is also reflected in the pseudoscalar one, in such a way that the difference between the 
two remains smooth. We have evaluated this difference within our model and it is also shown in fig. 3. It starts to 
decrease linearly with the density but then flattens off, due to the action of the nucleonic reaction which limits the 
rise with density of the mean scalar field. The convergence effect between the two susceptibilities is mild : with both 
quantities normalized to the vacuum scalar one the difference between the two susceptibilities is 28 in the vacuum 
and it becomes 23 at po- This evolution reflects the restoration associated with the smooth drop of the condensate. 
However the magnitude of the convergence effect is less marked that previously obtained in our previous work "^l. 
This difference with our previous conclusion is due to the interpretation of the condensate in the expression of the 
pseudoscalar susceptibility. The scalar susceptibility term of the quark condensate was not included in our previous 
work. 



IV. CONCLUSION 



In summary we have formulated a relativistic effective theory of nuclear matter based on the standard chiral 
symmetry properties. We have integrated in our approach the concept of the nucleonic response to the nuclear scalar 
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field, which is necessary to reach a reasonable description of the saturation properties of nuclear matter, working 
at the Hartree level. Our present mean-field description ignores the role of the pion loop. For a more complete 
description the corresponding correlation energy has to be included and we expect some modification on our values 
of the parameters such as the bare sigma mass. 

In our framework we have also studied the relativistic RPA, scalar or vector isoscalar, responses of the nuclear 
medium. The residual interaction is strongly density dependent. It turns from attraction into repulsion in the vicinity 
of the normal density, hence changing the collective nature of the responses. The coupling of the quark density 
fluctuations to the nucleonic ones makes the QCD scalar susceptibility sensitive to the collective behavior of the 
nuclear response. It is appreciably enhanced in the low density region where the residual ph force is attractive. This 
collective behavior of the nucleon-hole polarization propagator is also reflected in the pseudoscalar QCD susceptibility. 



APPENDIX A: RELATIVISTIC BARE POLARIZATION PROPAGATORS 



First consider the bare first order polarization propagator in the vector channel : 

nyy(O) _ f Sd^P QiPF - p) Qi\p + q\ - Pf) 



= - lim / ' (Al) 



We multiply the denominator and the numerator by Et^^ + E* whch leads to : 



Kvi^) _ v.^ f »d^p ^ ^,,eiPF-p)ei\p + q\-PF) 



hm / — ^ {Et^^ + E;) \ (A2) 

gi Q^oj {2TTy ^^'^ ^ {p + q )-p 

Obviously, the factor {E^^^ + E*) in the integrand can be replaced by 2 Ep since in the limit q^Q, p has to lie on 
the Fermi surface. It follows that : 

n°vv(0) ^ _E]^ f 8fp QipF~p)e{\p + q\-pF) . 



gl M*Ng^oJ {2nf {p + f)/2M*^ - p^/2M 



N 



and the remaining integral is manifestly the non relativistic no(0). This establishes the second relation given in eq. 
1221 For each scalar vertex one gets an additional multiplying factor M^/Ep. Thus, for one scalar vertex as is the 
case in the mixed term ngy(O), the relativistic correction factor altogheter disappears, hence the relation in eq. |221 
The presence of two scalar vertices leads to the first relation for the pure scalar polarization propagator n|5.g(0). 
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